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Abstract 

Random diffeomorphisms with bounded absolutely continuous noise are known to possess 
a finite number of stationary measures. We discuss dependence of stationary measures on an 
auxiliary parameter, thus describing bifurcations of families of random diffeomorphisms. A 
bifurcation theory is developed under mild regularity assumptions on the diffeomorphisms and 
the noise distribution (e.g. smooth diffeomorphisms with uniformly distributed additive noise 
are included). We distinguish bifurcations where the density function of a stationary measure 
varies discontinuously or where the support of a stationary measure varies discontinuously. 

We establish that generic random diffeomorphisms are stable. Densities of stable stationary 
measures are shown to be smooth and to depend smoothly on an auxiliary parameter, except at 
bifurcation values. The bifurcation theory explains the occurrence of transients and intermit- 
tency as the main bifurcation phenomena in random diffeomorphisms. Quantitative descriptions 
by means of average escape times from sets as functions of the parameter are provided. Further 
quantitative properties are described through the speed of decay of correlations as function of 
the parameter. 

Random endomorphisms are studied in one dimension; we show that stable one dimensional 
random endomorphisms occur open and dense and that in one parameter families bifurcations 
are typically isolated. We classify codimension one bifurcations for one dimensional random 
endomorphisms; we distinguish three possible kinds, the random saddle node, the random ho- 
moclinic and the random boundary bifurcation. The theory is illustrated on families of random 
circle diffeomorphisms and random unimodal maps. 
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1 Introduction 

To fix thoughts, consider a single map with an attracting fixed point P. Write W^{P) for its basin 
of attraction. Adding uniform noise of small amplitude gives a random map with a stationary 
density with support near P. Increasing the amplitude of the noise leads to a bifurcation when 
orbits can escape from W^{P). Two possibilities occur: escaping orbits can or cannot return near 
P. Escaping orbits lead to transient dynamics if orbits do not return or intermittent dynamics 
if orbits do return. How can such transitions occur? What are the quantitative characteristics? 
As a second issue, take a one parameter family of maps that exhibits bifurcations and add small 
bounded noise to it. What happens to the bifurcation set and in what way do bifurcations in the 
randomly perturbed map manifest? More generally, one can consider maps and families of maps 
where noise is an intrinsic part of the description. 

It is the purpose of this paper to work out bifurcation theory of random smooth diffeomorphisms 
from the perspective of stationary measures, providing answers to questions like the ones just stated. 
The context in which we perform this study is that of points being mapped into bounded domains 
according to a probability distribution, under some regularity conditions. One can think of points 
being mapped by a diffeomorphism, defined on some compact manifold Ai, followed by a random 
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perturbation. This defines a discrete Markov process on M given by transition functions P{x,A) 
providing the chance that a point x ^ M. ends up in a Borel set A C We will assume that the 
region where x is mapped into, is a bounded domain Ux- We argue that from a modeling point of 
view there are clear and good reasons to consider bounded noise: in most physical systems random 
perturbations are limited in their effect. 

We indicate the regularity conditions assumed in this paper. A first regularity assumption is 
smoothness of the density of the transition functions P{x,A). It is not assumed that this density 
vanishes on the boundary of its support; densities that are positive on bounded domains to model 
uniform noise are incorporated. For y G Ad, write Vy for the set of points x G Ad that are mapped 
to domains that include y; 

Fy = {x G Ad I y G support P{x, •)}. 

As a second regularity assumption, we suppose that the sets Vy are domains with piecewise smooth 
boundary varying smooth with y. This assumption is natural in the context of diffeomorphisms 
followed by noise, but does not hold in the context of endomorphisms (possessing critical points) 
followed by noise. Precise formulations follow in Sectiou fl . 1 1 below. 

An alternative description of the setup is by starting with a collection T of maps on Ad and 
a measure on T. Maps are drawn randomly, and independently, from T according to the given 
measure. Similarly one can consider maps depending on parameters that are drawn randomly. The 
random parameters are drawn from a bounded domain according to a given distribution. Typical 
examples are given by smooth maps / with additive or parametric noise. It will turn out that there 
is no loss of generality, as far as statistical descriptions are concerned, when considering maps with 
finitely many random parameters. We address an appendix to the exploration of the range and 
connections of these definitions. 

In the following we will mostly speak of random maps; maps depending on finitely many param¬ 
eters that are drawn randomly. In the part of the paper developing the general theory. Sections [3 
toEl we assume the maps to be diffeomorphisms. This guarantees the regularity assumptions for¬ 
mulated above in the description as discrete Markov processes to hold. Random endomorphisms 
are studied in one dimension in Section |H1 

Under the mild regularity assumptions, random diffeomorphisms possess finitely many station¬ 
ary measures, whose support is the closure of an open set and whose density functions (stationary 
densities) are smooth on all of Ad. The stationary densities are flat along the boundary of their 
support. This has immediate consequences for the statistical properties of orbits: orbits are very 
rarely found near the boundary of the support of the stationary densities. The main focus of this 
paper lies then in the description of the dependence of stationary densities on the random diffeo¬ 
morphisms. This includes describing quantitative characteristics. Motivated by examples discussed 
below we call a random diffeomorphism stable if its stationary densities and their supports vary 
continuously with the random diffeomorphism (precise definitions follow shortly). Otherwise we 
speak of a bifurcation. 

Below we introduce the precise setup and present our main results in a series of theorems. The 
presentation of the material is separated in a section treating random diffeomorphisms fSection ll .1 j) 
and a section treating families of random diffeomorphisms fSection 11 . 21 ) . Various aspects of the 
bifurcation theory for random diffeomorphisms are presented in Theorems II ..SI II .til II .71II .1 0111 .LSI 
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and ll.ldl Appendix comments on the setup. The main body of the theory is developed in 
Sections mtoci Sections mEi and |1] contain respectively material on transfer operators, proofs of 
stability theorems, and discussions of parameter dependence. Sections [HI and IHl develop material 
on conditionally stationary measures and apply this to compute expected escape times. Section [7| 
treats the speed of decay of correlations depending on a parameter. Appendix^contains an implicit 
function theorem used to obtain regularity in the parameter of solutions of integral equations 
involving the transfer operator. 

There is no easy analogous theory for random endomorphisms; stationary measures for random 
endomorphisms will in general be less regular resulting in different statistical properties. In Sec¬ 
tion El we present a satisfactory theory for one dimensional random endomorphisms, including a 
classihcation of possible bifurcations. This extends the general theory in particular by classifying 
codimension one bifurcations. Typically only finitely many bifurcations occur in one parameter 
families of random endomorphisms, in contrast to families of deterministic maps. The material on 
one dimensional random endomorphisms is developed and presented in Section El 

Section M contains two worked out examples of random circle diffeomorphisms and random 
unimodal maps. In Section fh.ll we consider the standard circle diffeomorphism with small additive 
noise; 

fa{x;uj) = X + a + Lo + esm{2TTx) (1) 

with X on the circle M/Z, for fixed e G (0,1) and a random parameter uj taken uniformly from a 
small interval. Section discusses as a prototypical example of random endomorphisms on an 
interval, random logistic maps 

f{x;i^) = {n + u})x{l-x) (2) 

with multiplicative noise obtained by varying uj with a uniform distribution in some interval. 

There is a large body of literature on stochastic stability (e.g. [SSIEZIIZIIIO], see also HU), 
considering bounded noise as a means to treat properties of single deterministic systems. This is 
done by letting the noise level decrease to zero. In contrast, we consider maps and families of maps 
where noise is an intrinsic part of the description. 

Previous attempts to study stochastic bifurcations fall into two categories. One way is to 
consider notions close to the traditional understanding of bifurcations for deterministic dynamics 
by embedding the random dynamical system into a skew product system. Another way studied 
frequently in the literature uses singularity theory to describe changes in the density of a stationary 
measure. This approach has been used to study systems with unbounded noise, often of a Gaussian 
nature, so that a unique smooth stationary density occurs. See p] for further discussion. We 
consider the shape of the stationary densities, but the noise being bounded leaves dynamics part 
of the picture. 

Central in control theory are control sets, that is maximal sets of approximate controllability. 
The shape of the control sets varies with external parameters, where discontinuous changes are 
possible. Interpreting the control variable as noise, a relation with the present paper becomes 
apparent. Discontinuous changes of the control set, interpreted this way, are among the bifurcations 
identified in this paper. In the context of differential equations depending on a control variable 
the study of control sets is taken up in [mEii. Bifurcation theory for such random differential 
equations in the spirit of this paper is considered in m- 
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1.1 Random diffeomorphisms 

The adjective smooth stands for Let be a smooth n-dimensional compact Riemannian 

manifold with measure m induced by the Riemannian structure. Let A be a closed domain in 
n-dimensional Euclidean space. Smoothness of a function <7 on A is to be understood in the sense 
that g can be extended to a smooth function on a neighborhood of A. 

Definition 1.1 A smooth random map, or random endomorphism, is a smooth map / : At x A —> 
AA, X I—> /(x;a;), depending on a random parameter lv G A drawn from a measure on A with smooth 
density funetion o' : A —> M, w 1 —> g{uj). A random diffeomorphism is a smooth random map so that 
X I—> f{x]io) is a diffeomorphism for each to. 

Remark 1.2 Alternatively one can explicitly include the noise distribution g as part of the defini¬ 
tion of smooth random map (and speak of a pair {f,g)). For convenience we consider g given. The 
results in this paper have direct, easily obtained, analogs if g is allowed to vary. 

Note that endomorphisms and diffeomorphisms are always assumed to be smooth. The basic 
setup we are treating is of points being mapped into bounded domains according to some probability. 
The following standing assumptions will be made with this setup in mind. The random parameters 
will be chosen from a region A that is a domain in M” with a piecewise smooth boundary. The 
number of random parameters is in particular equal to the dimension of the state space AA. The 
most important examples are where A is the unit ball A = {x G M” | ||x|| < 1} or the unit box 
A = {(xi,... ,Xn) G M” I |xi|,..., \xn\ < !}■ Throughout this paper we assume that to 1 — f{x]u>) 
is an injective map for each x. Hence /(x; A) is diffeomorphic to A. 

Write o for the measure on A with density function g. A smooth random map gives rise to a 
discrete Markov process through the transition functions 

P{x,A)= f dv{uj) (3) 

J { 1 ^ 1 f(x;uj)eA} 

for Borel sets A. With hx{io) = /(x;a;), the measure P{x, •) equals {hx)*n defined by {hx)*r'{A) = 
i'{hf^{A)). Vice versa, a discrete Markov process with noise from a ball or a box such that its 
transition functions have smooth positive densities admits a representation by smooth random 
maps (depending injectively on a random parameter), see Appendix El Some cases of parametric 
noise, where the maps do not depend injectively on the random parameter, do not fall into this 
setup. We refer to AppendixEIfor further discussion. 

The general theory will be developed for random diffeomorphisms, instead of random endo¬ 
morphisms. Endomorphisms allow for pathological examples, for instance maps f{x',oj) that are 
constant in x. We will however discuss random endomorphisms in one dimension (on a circle or a 
compact interval) in detail. 

With a slight abuse of notation, iterates of /(x;a;) are given as 

/^(x;wi,.. .,u}k) = f{f’'~^{x;iOi,... ,uJk-i)-,uJk). (4) 
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More generally, write for all infinite sequences u = {uji}i>i with each Wj G A. Denote f^{x; u) = 
... , 0 ;^). Let d : A^ ^ A^ be the left shift operator; = {uJi}i> 2 . Consider the 

skew product system S \ M. x A^ A4 X A^ given by 

S{x,u:) = {f{x-uji),'duj). (5) 

On A^ one considers a measure which is the product of the measure u over each A. 

With these definitions in mind, we introduce the central notions of stationary measures and 
ergodic measures. A stationary measure /x for the smooth random map / is a probability measure 
on Ad with ^ x S-invariant; 


/i X u°^{S-\B)) =/i X 

for Borel sets B C M x A^. Equivalently, see EZHU, 

/i(A) = j P{x,A)dfj,{x) 

Jm 

for Borel sets A C Ad. We refer to a stationary density as the density of an absolutely continuous 
stationary measure. 

A stationary measure ^ is called ergodic if /i x 12 °° is an ergodic measure for S in the usual 
sense that invariant subsets of Ad x A^ for S have zero or full measure. See 123 for equivalent 
formulations. The Birkhoff ergodic theorem tells that for an ergodic stationary measure, 

Hf{x]u;))= [ (j){x)dfj.{x) (6) 

JM 

for all integrable functions (/> on Ad and /r x almost every point {x,u). Taking (j) = 1 a, the 
characteristic function of a Borel set A C Ad, it shows that the relative frequency with which 
typical orbits visit A is given by /r(A). 

Write i?^(Ad) for the space of random diffeomorphisms / on Ad (with f{x',Lo) jointly in 
X € Ad and cj G A), depending on a random parameter from A through a distribution with a 
density function g. 

Let a random diffeomorphism / G i?°°(Ad) be given. The existence of finitely many ergodic 
stationary measures for / presented in the following theorem, can be found in HU Chapter 5] (valid 
under more general conditions). Similar results are contained in [1]. We add statements on the 
regularity of the stationary measures valid in our context. Differentiability of stationary densities 
is also discussed in imiHi. 

Theorem 1.3 The random diffeomorphism f G i?°°(Ad) possesses a finite number of ergodic sta¬ 
tionary measures pi,..., pm with mutually disjoint supports Ei,... ,Em- All stationary measures 
are linear combinations of pi,, pm- 

The support Ei of pi consists of the closure of a finite number of eonneeted open sets Cl,... ,Cf 
that are moved cyclieally by /(•; A). The density fii of pi is a C°° funetion on Ad. 


lim — 

n— >00 n 

1=0 
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Proof. See Chapter 5] for the existence proof of cyclically permuted ergodic stationary 
measures. 

Since a point x is mapped to a set /(x;A) diffeomorphic to A, the support Ei is the closure 
of finitely many connected open sets. We claim that the closures of the sets Ei, 1 < i < m, are 
mutually disjoint. Suppose on the contrary that dE^ and dEj with j ^ i have a point z in common. 
Then 2: is mapped by /(•; A) to the injective image of A. By invariance of Ei and Ej under /(•; A), 
f{z', A) is contained in both Ei and Ej, which is not possible. 

The regularity statements follow from Proposition 12.111 in Section |2 □ 

Densities of stationary measures for random endomorphisms are in general not smooth functions, 
but are less regular (see Section El- Random endomorphisms without critical points, such as 
expanding maps, do possess finitely many smooth stationary densities. The proofs for random 
diffeomorphisms extend to cover such random endomorphisms. The regularity of a stationary 
density (j) implies that cf) is flat along the boundary of its support E. By the Birkhoff ergodic 
theorem (applied to the characteristic function of a neighborhood of BE), this means that typical 
orbits are very infrequently found near the boundary of E. 

We introduce a topology on the space R^{A4) of random diffeomorphisms in order to be able to 
compare the dynamics of nearby random diffeomorphisms. Natural topologies on R^{A4) are the 
uniform topologies on C’^iM X A,M). See e.g. EH for generalities on these topologies. We 
will assume R^{A4) to be equipped with this topology. Note that the alternative approach through 
discrete Markov processes suggests a topology using the densities of the transition functions. 

Consider / € R°°(A4). Write /xi,... ,/Jm for the stationary measures of / G R°°{M.) given by 
Theorem o 

Definition 1.4 R random endomorphism f € R°°(j\4) is stable if for all f sufficiently close to f, 
the following two properties are satisfied. 

• For each i,l < i < m, the random endomorphism f has a stationary measure fii whose 
density is close to that of pLi. 

• The supports of fii and p,i are close in the Hausdorff metric. 

We speak of a bifurcation, or a bifurcating random endomorphism, if at least one of these properties 
is violated. 

Definition 1.5 An ergodic stationary measure of f (z R°°{M) is called isolated or attracting, if 
there exists an open set W (an isolating neighborhood) containing the support E of pL, so that 
f{W;A) C W and pL is the only ergodic stationary measure of f with support in W. 

For each / close to / is the closure of f{W;A) contained in W. The following stability result 
shows that nearby random diffeomorphisms have indeed a unique stationary measure with support 
in W. For bifurcations where stationary measures vary discontinuously the condition of being 
isolated must therefore be violated. The proof of the following theorem is found in Sectional 
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Theorem 1.6 Let /x be an isolated ergodic stationary measure of f ^ R°°{M.) with density 4> with 
isolating neighborhood W. Then each f G sufficiently close to f possesses a unique ergodic 

stationary measure ft with support in W. The density (j) of fi is C°° close to f. 

Note though that the above theorem leaves open the possibility that the supports E oi fi and E 
of /i are not close in the Hausdorff metric. An illustrative example of this phenomenon is described 
in Section EH 

The following theorem establishes that stable random diffeomorphisms are generic. Its proof 
is in Section E The argument also shows that random diffeomorphisms with a locally constant 
number of smoothly varying stationary densities (ignoring variations in their support) form an 
open and dense subset of R^{Ai). 

Theorem 1.7 The set of stable random diffeomorphisms in R°°[M) contains a countable inter¬ 
section of open and dense sets. 

A thorough description of the dynamics of random circle diffeomorphisms is in Section Ih.l I 
Random diffeomorphisms, and even random endomorphisms, on the circle are shown to form an 
open and dense set, see Theorem 18.hi 

1.2 Families of random diffeomorphisms 

Bifurcations are best studied in families depending on finitely many parameters. We will consider 
families of random diffeomorphisms depending on a single real parameter, where we have the goal 
to focus on bifurcations that typically occur varying one parameter. 

Definition 1.8 A smooth family of random endomorphisms is a family of random endomorphisms 
{fa} depending on parameters a, so that fa{x-,Lj) depends smoothly on {x,uj,a). A smooth family 
of random diffeomorphisms is a smooth family of random endomorphisms where each map fa{-]co) 
is a diffeomorphism. 

Remark 1.9 Alternatively, one can explicitly include noise densities ga in the definition (consider¬ 
ing pairs {fa, da)) with ga{w) varying smoothly with {uj,a). Compare Remark \l.E. For convenience 
we consider fixed noise densities, but completely analogous results hold if the noise densities are 
allowed to vary with a. 

Consider a smooth one parameter family {fa} of random diffeomorphisms, with a from an 
interval I. Consider a parameter value oq and an ergodic stationary measure HaQ with support 
EaQ. The following result extends Theorem II .61 providing an analogous statement in the context of 
families. If pao is an isolated ergodic stationary measure then there are ergodic invariant measures 
Pa for o near oq with nearby densities. 

Theorem 1.10 Suppose pao is an isolated ergodic stationary measure. Then the stationary density 
X I—> (fa{x) of Pa depends C°° on {x,a). 


See Section 0 for the proof. We stress again that the support Ea of fia can still vary discontinu- 
ously in the Hausdorff metric with a. The number of components of the support of the stationary 
measure can also change, while the stationary density varies smoothly. 

Consider a smooth function cj) with support on an isolating neighborhood W for fiao and compute 
averages of (j) along orbits /„ (x; u). By the Birkhoff ergodic theorem, for typical initial points x 
and noise sequences u, the averages lie on a smooth function of a for a near oq. 

The two types of bifurcation distinguished in Definition 11.41 gives rise to a particular dynamical 
phenomenon associated to either intermittency or transients. Consider a family {fa} of random 
diffeomorphisms in R°°{A4), with a from an open interval I. Suppose that oq G / is a bifurcation 
value for {fa} involving a stationary measure //. Write (p for the density of /i. Analogies with 
deterministic dynamics suggest the following two definitions. In Section |H1 we will see that in 
typical one parameter families of random interval or circle endomorphisms bifurcations are isolated 
and of these two types. 

Definition 1.11 The bifurcation at oq is called an intermittency bifurcation if there is a stationary 
density (pa for {fa} with pao = P and depending continuously on a, so that the support Ea of pa 
varies with a, for a near oq, as follows. 

• Ea varies continuously for a from one side of oq . Without loss of generality, we assume this 
to be the case for a < oq. 

• Ea is discontinuous at a = oq and Ea contains an open set disjoint from Ea^ for a > oq. 

An orbit piece outside a small neighborhood W of Ea^ is called a burst. Out of the substantial 
literature on intermittency in dynamics, we point to references iiniiiiiisniiiHiisiiEsi- 

Definition 1.12 The bifurcation at oq is called a transient bifurcation if there is a stationary 
density pa for {fa} with pa^ = p for a close to oq from one side of oq (without loss of generality, 
we assume this to be the case for a < oq), so that 

• pa and its support Ea vary continuously with a, for a < oq. 

• there is no stationary density near pa^ for a close to oq and a> oq. 

We end this section with a quantitative description of time series near intermittent or transient 
random bifurcations. We do this through the estimation of the expected escape time from a small 
neighborhood of the support of the bifurcation stationary measure. Estimating the speed of decay 
of correlations as function of a parameter gives further details of changes through bifurcations. 

As before, let {fa} be a family of random diffeomorphisms in with a from an open 

interval I. Let fiao be a stationary measure of /„□ for some oq € I and let W be an open neigh¬ 
borhood of the support Eag of flag such that no other stationary measure has support intersecting 
W. If fag is stable, then there is a unique stationary measure fia that is the continuation of Hao- 
The stationary measure fia has its support in W for a near oq. If oo is a bifurcation value for {fa}, 
iterates fa{x\u>), for certain x G IE and a near oq, may leave IE. 
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Figure 1: Typical times series for intermittent dynamics on the left and transient dynamics on the right. 
The time series are for bifurcation values after the bifurcation took place. The intermittency bifurcation 
involves interval diffeomorphisms with a single stationary measure; the support consisting of two intervals 
for a period two cycle bifurcates to form a single interval. In the transient bifurcation one stationary measure 
out of two stationary measures existing previous to the bifurcation disappears. 



For X G PF and u G A^, define 

Xa{x,u) = mm{k \ /a(x;cj) ^ W). (7) 

The following result shows how the average escape time from a neighborhood of the support of a 
bifurcating stationary measure is more than polynomially large in an unfolding parameter. This 
makes it difficult to accurately establish the bifurcation parameter value using finite data, even in 
numerical simulations. It explains the occurrence of very long transients near a transient bifurcation 
and the very irregular occurrence of bursts in intermittent time series. The proof, in Section El 
relies on the construction of conditionally stationary measures in Section El 


Theorem 1.13 For each k > 0 there is a constant Ck > 0 so that 



Xa{x,co)di'°°{u)dm{x) > Ck\a — oq 


A single random map with an isolated measure supported on a single component has exponential 
decay of correlations as precised in the following proposition. The interest from our perspective in 
computing the speed of decay of correlations lies in the study of bifurcations where the support of 
a stationary measure has several components merging. This will be discussed below. Proofs of the 
following statements are in Section [3 The reader can consult mM for background on decay of 
correlations. Write 



'll; o p{x-,uji,... ,u;n)diy{u;i) ■ ■ ■ dv{ijjn)- 
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Proposition 1.14 Let f be a random map with an isolated stationary measure p. with connected 
support. Let W be an isolating neighborhood for fi. Take ^ G C?‘fW'). Then 



ip{x)U'^'il){x)dm{x) 



ip{x)dm{x) 



ijj{x)dp{x) 


< 


for some C > 0, 0 < p < 1. 


Note that the exponential decay of correlations holds for observables (/?, V' £ C^iW), thns in¬ 
cluding characteristic functions of open sets. On the other hand, our definition involves an average 
over noise sequences. The following remark addresses this point. 


Remark 1.15 For a fixed noise sequence u € and nonnegative observables '0 G C'^iW), 
consider 



ip{x)if o p{x;uii,... ,ujn)ix)dm{x) 



ip{x)dm{x) 



ijj{x)dfi{x) 


( 8 ) 


By Proposition \1.14[ the integral over A” of this expression is bounded by Cp'^. Necessarily, ^ 
is exponentially small in n for {uji,... ,u}n) outside an exponentially small set. Indeed, choose p 
slightly larger than p. Then is larger than p'^ only on a set Sn C A"' with u'^{Sn) < Cp'^/p'^. 


We return to a family {fa}, a G /, of random maps. Assume that fa has an isolated measure 
Pa for all a G / with an isolating neighborhood W. Suppose qq G is a bifurcation value for an 
intermittency bifurcation so that 

• the support of pa consists of k components for a < ag, 

• the support of pa consists of a single component for a > ag- 

We incorporate the dependence of t/” on a into the notation by writing [/”. 

Theorem 1.16 Let {/„} be as above. Take € C^iW). There are a constant C > 0 and a 
smooth function a ^ pa with Pa = ^ for a < ag and Pa < ^ for a > ag, so that 


for a > ag. 



ip{x)U 2 'f’{x)dm{x) 



ip{x)dm{x) 



if{x)dpa{x) 


< Cp 


n 

a 


The smoothness properties of pa imply that 77a — 1 is a flat function of a at a = oq. As in 
the discussion of escape times, that shows how slowly the bifurcation manifests itself in time series 
when moving the parameter a. Remark 1 1.15 1 also applies in the parameter dependent context. 
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2 Transfer operators 


Let / be a random diffeomorphism on the manifold A4. Associated to / is the stochastic transition 
function, 

P{x,A)= [ dv{uj) = [ lA{f{x]u;))di'{uj), 

J{lo I f{x-,uj)GA} Ja 

for Borel sets A C Ai. Write hx{uj) = f{x]uj), note that h maps A injectively onto Ux = f{x]A). 
The density 

k{x,y) = d{hx)^u/dm (9) 

of P{x,-) vanishes for y outside Ux and is a smooth function on its support. 

Write £^(A4) for the space of integrable functions on Ai. Dehne the transfer operator L acting 
on C^{Ai) by 

/ L(j){x)dm{x) = / P{x, A)(l){x)dm{x). (10) 

Ja Jm 

The transfer operator is a positive linear operator. A stationary density is a fixed point of L. 
Define 

Vx = {z € Ai \ X e f{z] A)}, (11) 

which is the set of points in Ai that are mapped to x by some random map. The assumptions 
on the random diffeomorphism / imply that the equation f{x\uj) = y can be solved for x as a 
diffeomorphic map of u). Therefore Vx is diffeomorphic to A and thus a domain with piecewise 
smooth boundary, depending smoothly on x. 

Write Pf(.-uj) for the Perron-Frobenius operator, defined by 


Pf(-,u)4>{x)dm{x) = 




(j){x)dm{x) 


( 12 ) 


for Borel sets A. That is, for a measure y, with density (p = dy/dm, f^y has density Pfp = df^^y/dm 
(see e.g. m)- The following lemma gives the transfer operator as an average over the random 
parameters co of the Perron-Frobenius operators for and gives an equivalent formulation as 

an integral over the state space Ai. 


Lemma 2.1 The transfer operator is given by 


L(/>(x) = / Pfy^)p{x)diy{Lj), 

JA 

(13) 

or. 


L(j){x)= k{y,x)(j){y)dm{y) 

(14) 

■Jv^ 


Proof. By (IIOI). for a continuous function ip, 




'ip{x)L(l){x)dm[x) 




ip{f{x', uj))(p{x)dv{ijj)dm{x) 
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Calculate 




'0(/(x; uj))(j){x)dv{u})dm{x) 


This implies m- Alternatively, 


/ / '>P{f{x]uj))(j>{x)dv{uj)dm{x) 

Jm Ja 


proves m- 



A JM 



A JM 


'tp{f{x‘, uj))4>{x)dm{x)di'{uj) 
’ip{x)Pf(^..^4>{x)dm{x)di'{uj) 


IM Ju, 


'M JVy 


'ip{y)(p{x)k{y, x)dm{y)dm{x) 
'ip{y)(p{x)k{y, x)dm{x)dm{y) 


□ 


Remark 2.2 The transfer operator L preserves integrals, as the following computation shows. 

/ Lcl){y)dm{y) = / / k{x,y)4>{x)dm{x)dm{y) 

Jm Jm Jvy 

/ k{x,y)(j){x)dm{y)dm{x) 

JUrr. 


IM JU,y 

/ (l){x)dm{x) 

iM 


because k{x,y)dm{x) = 1. 
Iterating L gives 

L^(^(x) = 


'14 


k{z, x)L(f){z)dz 


'V4 JV. 


k{z,x)k{y,z)4>{y)dydz 


lf-^V,,A) -114 


k{z,x)k{y,z)dz(j){y)dy 


which is of a similar form, namely a) x)4i{y)dy with k 2 {y,x) = Jy k{z, x)k{y, z)dz, as 

L(j){x). Inductively similar expressions are derived for higher iterates of L. 

Denote by V{M.) = {(f) ^ C^{M) | > 0, ^M 4>{x)dm{x) = 1} the space of densities on A4. The 

above remark shows that L maps V into itself. Smoothness of stationary densities is obtained by 
showing that L maps a space of smooth densities into itself. 


Proposition 2.3 The transfer operator L maps C^{A4) into itself and is a compact operator on 
C^(M). 

The number 1 is an eigenvalue of L with equal algebraic and geometric multiplicity m > 1. 
The densities (pi,... ,(j)rn provide a basis of eigenfunctions with mutually disjoint support. Each 
eigenfunction (pi is C°° and its support consists of a finite number Ci of connected components. 
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Proof. Theorem gives the m invariant densities (pi,..., 4>m- The geometric multiplicity of 
the eigenvalues 1 is equal to m. Since L preserves the norm, the algebraic and geometric 
multiplicity of 1 are equal. To see this, suppose on the contrary that there is a nontrivial vector 
in ker(L — /)^\ker(L — I). Elementary linear algebra gives the existence of a sequence of vectors 
^pn inside ker(L — I)^ converging to an eigenvector (p, such that lim^^oo = ‘^(p- Indeed, take 

pj € ker(L — I)^ with Lip = cp + ip and let ipn = ^ip- From L'^{ip) = n(p + ip it follows that 
L'^{(p + Ipn) = 2(p + Ipn, which converges to 20 if n —> oo. This contradicts the preservation of the 
norm by L. 

There can be no additional eigenvectors of L that do not correspond to linear combinations of 
densities. Namely, suppose (p is an eigenvector taking both positive and negative values. Write 
(p = (p^ — (p- for nonnegative functions 0+ = max{0,0}, 0_ = max{—0,0}. If 0 is not a linear 
combinations of densities, the supports of 0+,0_ cannot be invariant. Since L is positive and 
preserves the norm, (T0)_|_ < T0_|_ so that 0 cannot be an eigenvector. 

We will show that L maps into C^(A4) and C^{A4) into From this it follows 

that (pi € Take ip € Use a chart to identify a neighborhood of x € A4 with an 

open set in M”. With h a small vector in M”, consider 

Lip{x + h) — Lip{x) = 


The first term on the right hand side is small for h small by continuity of k and integrability of ip. 
The other two terms are small for h small by the continuous dependence of Vx on x. Continuity of 
Lip follows. Suppose next that ip G and consider {Lip{x + h) — Lip{x)) This equals the 

right hand side of divided by \h\. Note that 

+ h) - k{y,x))ip{y)dm{y) = f ■^k{y,x)-l^ip{y)dm{y) 

is a continuous function of x. To check continuity of the remaining two terms it suffices to do a 
local calculation by covering the boundary of 14 by finitely many balls and using a partition of 
unity. Without loss of generality we may assume that near a smooth part of the boundary of 14) 
14 is bounded from below by the graph of a continuously differentiable function Hx ■ [0,1]”“^ —> M. 
We may also assume that m equals Lebesgue measure on M”. Then 


/ k{y,x)ip{y)dm{y)- k{y,x)ip{y)dm{y) 
/ {k{y,x + h)-k{y,x))ip{y)dm{y) 

I 


+ / k{y,x + h)ip{y)dm{y) 


'U\(U+;,nl4) 


k{y,x)ip{y)dm{y). 


(15) 


h-'O \h\ 


1 r rHoo+hiyi) 

k{y,x)ip{y)dy = lirn — / / k{yi,y 2 ,x)ip{y)dy 2 dyi 

f d h 

= / ■^Hxiyi)-rT^k{yi,Hx{yi),x)ip{yi,Hx{yi))dyi 
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is a continuous function of x. The contribution near the finitely many points where I4 is not 
smooth vanishes in the limit /i —> 0. Summarizing, D[L'ip) has an expression of the form 

D{L'il;){x) = [ ■^k{y,x)Tp{y)dm{y) + [ n{y,x)k{y,x)^p{y)dS{y) (16) 

Jv. ox Jav. 

where n{y, x) measures the change of dVx in the direction of the unit normal vector to I4 and dS 
is the volume on dVx- We remark that the formula is a variant of the transport theorem 7.1.12 and 
the Gauss theorem 7.2.9 in pp. It follows that Lip is continuously differentiable if ip is continuous. 
Higher order derivatives are computed inductively. This gives that Lip £ for ip £ C^{Ai). 

We prove compactness on C^{M) by modifying the argument in [3J. Let be the unit 

sphere in C^{A4). Consider first A; = 0. By the Arzela-Ascoli theorem, compactness of L on C^{A4) 
follows from the following two properties (compare |f)4j 1. 

• for all X £ M, {\Lip{x)\ \ ip £ is bounded, 

• LF is equicontinuous. 

For Ip £ B^(A4), \Lip{x)\ < k{y, x)dm{y). This is a continuous function of x and hence bounded. 
This proves the first item. The above computations showing that Lip is continuously differentiable 
also show that \\D{Lip){x)\\ is uniformly bounded on B^(A4). This proves that LF is equicontinu¬ 
ous. Compactness in C^(M) follows similarly by noting that 

• for all X £ M, i < k, {\\D\Lip){x)\\ \ ip £ B^{M)} is bounded, 

• \\D^~^^{Lip){x)\\ is uniformly bounded on il^(Ad). 

□ 

Remark 2.4 The transfer operator L is compact on the space of quadratic integrable func¬ 

tions on M. I.lill Section X.2] and The proof of Provosition \2JA demonstrating that the transfer 
operator increases regularity of functions, implies that the spectrum of L on equals that of 

L on C^{M). 

Remark 2.5 The spectral radius of L is 1, the eigenvalue 1 occurs with multiplicity m equal to the 
number of stationary measures. The peripheral spectrum on the unit circle consists of eigenvalues 
^ml-p ^ 0 < i for each p occurring as the number of connected components of a stationary 

measure. See and Theorem V.f.O] for a proof. 

Proposition 2.6 The transfer operator L as a linear map on C^{A4) or £^(A4) depends contin¬ 
uously on f & R^{M.). 

Proof. Consider / near /. Write L and L for the corresponding transfer operators. We need to 
prove that L—L has small norm. Consider the transfer operators operating on C^(A4) (continuity on 
£^(A4) is treated analogously). The transfer operator L is given as L(p{x) = jy k{y,x)(p{y)dm{y). 
For (p £ B^{Ai), the unit sphere in C^(A4), 

L(p{x) - L(p{x) = / k{y,x)(p{y)dm{y) - / k{y,x)(p{y)dm{y) 

Jv,, 
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is small, uniformly in x, since k{y,x) is close to k{y,x) on 14 n 14 and 14 is close to 14- The 
derivative D{L(j)) is given by (lltil) . An analogous formula holds for D{L(f)). Since the functions and 
sets involved in the two formulas for D{L(f)) and D{L4>) are close, D{L(j))(x) is uniformly close to 
D{L(j)){x). Closeness of higher order derivatives, up to order k, is treated analogously. Continuity 
on C‘^{AA) is proved analogously. □ 

Remark 2.7 Consider two nearby random diffeomorphisms f and f from R^{A4). Write L and 
L for the corresponding transfer operators on C^(A4). Let Xi,..., Xi be a finite set of eigenvalues 
for L and denote by F the sum of the corresponding generalized eigenspaces. Then L possesses a 
nearby set of eigenvalues Ai,..., A;. The sum F of the corresponding generalized eigenspaces is a 
small perturbation of F (in the sense that F and F have nearby bases). See l,V.5L Theorem IV.3.16]. 

3 Stable random diffeomorphisms 

This section contains the proofs of Theorem 11.61 on stability of isolated stationary measures and 
Theorem II .71 establishing generic stability of random diffeomorphisms. 

For the restriction of random maps to an isolating neighborhood W we consider the transfer 
operator acting on functions vanishing outside W and at the boundary of W. Write 

C'q(IT) = {/ G C^{M.) I the support of / is contained in W}. 

Then L acting on Cq{W) is well dehned. The results in the previous section hold for L acting on 
C^{W). 

Proof of Theor,em I1.6L Recall that the closure of /(VF; A) is contained in the isolating neigh¬ 
borhood W. This property extends to random diffeomorphisms sufficiently close to /. Restrict the 
map X I—> f{x]uj) to W and and consider the transfer operator L acting on C^fW). Then L has a 
single eigenvalue 1. Since the spectrum of the transfer operator varies continuously with the random 
diffeomorphism at /, the transfer operator corresponding to each nearby random diffeomorphism 
possesses a single eigenvalue 1. The corresponding eigenvector is near (]. □ 

Lemma 3.1 Write Lf for the transfer operator on C^(A4) for f G R°°{A4). Densities of station¬ 
ary measures vary continuously with f G R°°{M.) at a random diffeomorphism f precisely if the 
multiplicity of the eigenvalue 1 for Lj is locally constant in f for f near f. 

Proof. Consider / G R°°{M) with an eigenvalue 1 of multiplicity m. Let p,i,...,firn be the 
ergodic stationary measures with densities (j)i,... ,4>m- Write Fj be the direct sum of the lines 
spanned by (pi ,..., (pm. By Remark 12.71 the transfer operator for any / G R°°{A4) sufficiently close 
to / possesses a m-dimensional invariant linear space Ff that is the continuation of Fj. 

The spectrum oi Lf restricted to Ff is in general close to 1. Suppose now that all eigenval¬ 
ues equal 1. Then (p ^ Ff implies LfCp = (p. Write (p = (p~^ — (p~ with (p~^ = max{0, (/>} and 
(p~ = max{0, —(p} the positive and negative parts of (p. Because Lf > 0 and Lf preserves the 
norm, (Lf<pff < LfCp'^ precisely if (p~ and (p'^ are not invariant. Thus and p~ are necessar¬ 
ily invariant. It follows that invariant densities are obtained by taking positive parts of invariant 
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eigenfunctions. This way m invariant densities for / near those of / can be obtained, proving the 
lemma. □ 

Proof of Theorem II. 7L Consider diffeomorphisms on an open neighborhood 17 of a random 
diffeomorphism / € Write m for the multiplicity of the eigenvalue 1 for the transfer 

operator corresponding to /. There is a neighborhood ZD of 1 in the complex plane, so that for 
U small enough, each f ^ U has m eigenvalues counting multiplicity in D. Let F denote the 
m dimensional invariant linear space corresponding to these eigenvalues. Consider the map that 
assigns to f £ U the union of the support of all functions in F. By the continuous dependence of F 
on /, this is a lower semicontinuous set valued mapping and therefore continuous on a set B2 <Z U 
of Baire second category [24] . 

Consider the map that assigns to random diffeomorphisms / G the multiplicity m{f) of 

the eigenvalue 1 for the corresponding transfer map. By the continuous dependence of eigenvalues 
of the transfer map on /, the map m is upper semicontinuous. Since m takes on hnitely many 
values, it is continuous on an open and dense subset of R°°{M.). Indeed, consider An = {f £ 
R°°{A4) I m(/) < n}. The set of points of continuity of m, in the vicinity of some map in R°°{Jt4), 
equals the intersection of a finite collection of open and dense sets An U (Z?°°(ZW)\^„,), namely 
with n ranging over a finite set of positive integers. 

With reference to Lemma EH the two above items combined prove the theorem. □ 


4 Auxiliary parameters 

Consider a smooth one parameter family of random diffeomorphisms x 1 —> fa{x]U)) depending on a 
from an open interval I in M. The transition map P and its density k depend on a, we write Pa 
and ka- The support Ux,a of ka is assumed to vary smoothly with x and a. The density ka{x,y) 
is a smooth function of {a,x,y) £ Ua,x{a} x x Ux^a in the sense that it can be extended to a 
smooth function on an open neighborhood. Let La denote the transfer operator for fa, given by 

La(l){x)= kaiy,x)4>{y)dm{y). (17) 

The domain of integration I4,a = {u & Ai \ y £ fa{x', A)} depends smoothly on (x,a). 

Proposition 4.1 For r > 0, the transfer operator (j) 1 —> Laf) as a map from into 

is a Cr+i 

map of a and (f. 

Proof. By Proposition [THl La depends continuously on a. For the derivative of La/f with respect 
to a we hnd an expression similar to cni), 

d f d f 

—La(j){x)= —ka{y,x)(l){y)dm{y)+ Sa{y,x)ka{y,x)(j){y)dS{y) (18) 

da da Jov,,, 

for some smooth function Sa- It follows that for cf £ C^{A4), -^Laf) £ C^{Ai). This implies 
differentiability of (</>, a) 1 —> Lafj) for cf) £ Higher differentiability is treated similarly. □ 
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The operator La : C^{M.) —* C^{M.) does not depend on a, since -^La4> may not exist if 
(j) G C^{M.) and -^La4> is a function if </> G C^{A4). What does hold is that (x, a) i-H’ La4>aix) 
is if {x,a) 1 -^ 4>a{x) is . 


Proof of Theorem II.lOL Let W be the isolating neighborhood for ^ao- For a near uq, fa{W-, A) 
is strictly contained in W and fa has a unique stationary measure with support in W. Restrict fa 
to W for such values of oq. 

Consider the transfer operator La for /„ acting on Cq{W). Write F for the line in C^{W) 
spanned by 4>ao- Then Cq(W) = F®Hq{W) with Hq(W) consisting of functions with vanishing 
integral; 

H^{W) = {(/. G C^{W) I [ (P{x)dm{x) = 0}. 

Jm 


Write (fa for the eigenvectors of La continuing (fao provided by Proposition 11.61 Decompose (fa = 
if a with if a G Hq{W). Then if a is a solution of Laifa = a F (fao ~ La fan- Note that 




■ 


{x,a) ^ (faoix) -Lafaoix) is C"" 
Appendix I bI implies the result. 


The spectrum of L, 


ao 


Hf{W) 


is away from 1. Proposition ESI i 


m 


□ 


5 Conditionally stationary measures 

To study average escape times from open sets we make use of conditionally stationary measures, 
which are measures for which on average a fixed percentage of mass escapes under an iterate. We 
recall the notion of conditionally invariant measure, see gsi im cn ESI for its use in deterministic 
dynamics. Let a map f : A4 ^ A4 he given and restrict f to a domain W d M. Let V dW he 
the set of points in W that are mapped into W, points in the complement of P in kP are mapped 
outside W. Consider f '■ V ^ W. A conditionally invariant measure for / on VP is a measure on 
M so that n{A) = ^{f~^{A))/n{f~^iW)) for Borel sets A d W. 

Definition 5.1 Let f G Let W he an open domain in M. A measure jl on W is a 

conditionally stationary measure if 

fi{A) = f P{x,A)dfi{x) I f P{x,W)dfL{x) 

Jw / Jw 

for Borel sets A d W. 

See jlHl uni where this notion is called a quasistationary measure. Note that a conditionally 
stationary measure is a stationary measure if P{x, W)dft{x) = 1, that is, if the support of the 
conditionally stationary measure lies inside W. 

Lemma 5.2 A measure p, on W is a conditionally stationary measure for f if and only if p x v°° 
is a conditionally invariant measure for (/,■!?) on VP x A^. 

Proof. Write D^{x) = {a; G A^ | f{x;u>) G kP}. Consider S : \Jx^w{{x} x D^{x)) ^ kP x A^, 
S{x,Lj) = {f{x]Lo),Pu). We must show that the following two statements are equivalent. 
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(i) J1 X /ji X u°^{S-^{W X A^)) = /2 X u°°{A) for Borel sets ^ C x A^. 

(*0 Iw /a ^u{f{x-,uj))dfl{x)du{uj)/ lw(f(x;uj))d^(x)diy(cu) = lu(x)d^(x) for Borel sets 

U CW. 

Take a Borel set U x V with U cW and V C A^ and compute 
flxu°°{S-\U xV)) = 


Further 

fi X u°°{U xV) = [ lu{x)dfl{x)i^°°{V) (20) 

Jw 

Equations (dl) and (PH) contain the implication (i) =^> {ii) when applied for U x A^ for the 
enumerator and for W x A^ for the denominator. 

To show that {ii) implies (i), note that (fm and PH) show that {i) holds for Borel sets A of the 
form U X V if {ii) is assumed. Therefore it holds for all Borel sets in VF x A^. □ 

We continue with the introduction of transfer operators whose hxed points are the densities of 
conditionally stationary measures. The transfer operator Z, defined for functions in C^iW) with 
integral 1, is given by 

-^(</>) = / / L4>{x)dm{x). ( 21 ) 

/ Jw 

Write L(j) = 

Proposition 5.3 Z maps C^{W) into itself and is a compact operator on 
Proof. If k{x, y) denotes the density of the stochastic transition function 

L4>{x) = _ k{y,x)4>{y)dm{y). 

Jwnv^ 

Note that W nVx depends continuously on x. Recall from the proof of Proposition [231 that we 
must show that 

• for all X G VF, {|Z'0(x)| | G B^{W)} is bounded, 

• LF is equicontinuous. 


it. 

P{x, •), then 
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Here B^{W) is the unit sphere in C^{W). The first item follows as before: |L^/:(x)| < k{y, x)dm{y) 
is bounded by a continuous function and thus bounded. For the second item we must show that 
for each e > 0 there is 5 > 0 so that for all x € W, 'ijj € \L'ijj{x + h) — Lijj(x)\ < e if \h\ < 6. 

Recall, see lH^l) . 


L'iIj{x + h) - L'4){x) = / {k{y,x + h) - k{y,x))'ijj{y)dm{y) 

Jwnv^+hnV:, 

+ / + h)^jJ{y)dm{y) 

Jwnv^+h\{wnv^+hnV:c) 

- / Ky^x)'4){y)dm{y). 

Jwnv^\{wr\V^+hr\V^) 


Now 

1 / {k{y^x + h)-k{y,x))il){y)dm{y)\ < I \k(y, x + h) - k(y, x)\dm{y), 

Jwr\V^+hr\Vx 

which is small for \h\ small by uniform continuity of k. And 


1 / k{y,x+ h)'il;{y)dm{y)\ < j \k{y, x + h)\dm{y) 

Jwnv,+h\{wnv,+unv^) ./vunv;,+A(H^ni4+^nK) 

is small for \h\ small by boundedness of k and uniform continuity of the volume of 14 in x. Similarly 
for the third term. This proves equicontinuity. □ 


Recall from Proposition 12 .til that L depends continuously on the random diffeomorphism. The 
same argument shows that L depends continuously on the random diffeomorphism. 

Proposition 5.4 The transfer operator L as a linear map on C^{A4) depends continuously on 
f £R^{M). 

We obtain conditionally stationary measures by a perturbation argument, perturbing from 
an invariant measure. We do not develop general existence results for conditionally stationary 
measures, as such general results are not needed for our purposes. Let {fa} be a family of random 
diffeomorphisms depending on a real parameter a (z I. Consider, for oq £ I, a stationary density 
(pao with support Ea^ ■ Let VF be a neighborhood of Ea^ disjoint from the supports of possible other 
stationary densities of fa^- 

Proposition 5.5 For a close to oq, fa possesses a conditionally stationary density (pa on W, with 
4>ao = 0ao ond (x,a) i— (pa{x) continuous in (x,a). One has 

L(pa — (xlyO^Cpai 

where a(a) = La^a(x)dm(x) depends continuously on a. 
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Proof. The operator La varies continuously with a and therefore possesses a single eigenvalue 
close to 1 for a close to oq. The function (l)a is the corresponding eigenfunction. □ 


Recall the definition of the escape time Xa(a:,cj) forxGW and u G A^: 

Xaix,uj) = min{A; | /^(x;a;) ^ IT}. 

Lemma 5.6 

/ / Xa{x]OLi)di^°^{u;)(j)aix)dTn{x) = - - 

Jw Ja” 1 

Proof. Let Sa{x,uj) = (/a(x;ti;),'!?cj). Write 

D2 = {{X,UJ) I f2 {X-, UJi,... ,UJn) G W} 

for the set of points in VP x that remain in VP x A^ for n iterates of Sa- The exit set of 
points that leave VP x A^ in n iterates equals D^~^\D^. Thus Xa(x,u:) = n on E^. Write p,a for 
the conditionally stationary measure with density cj)a- From (11911 with A = VP x A^ we get 


a{a) = [ P{x,W)dnaix) 

Jw 


w JA 


lwif{x-, uo))dv{uj)d^Xa(x) 


= jxa^y°°{S-\Wxlx^)). 

It follows that fLa X — a). Calculate 

OO 

Xa{x-,uj)du°°{uj)dfla{x) = fla X (E^) 

k=l 

OO 



W JA^ 


k=l 

1 


— a 


1 — a 


□ 


As a corollary we obtain that the average escape time from VP goes to infinity as o —> oq. More 
precise estimates are derived in the following section. 

Proposition 5.7 Xa{x,^j^)diy°°{uj)dm{x) converges to oo as a oq. 

Proof. Let Ea be the interior of the support of (jia- 


Xa{x,oi))du°°{u)dm{x) > [ [ ^ dv°°{u)(i)a{x)dm{x) 

IwJA^ JEaJA’^ 4>a{x) 



> 


c [ [ Xa{x.,cj)dv'^{u)(j)a{x)dm{x) 

JEa JA^O 


c 


1 — a{a) ’ 
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for C = 1/max{x € W | 0a(a;)}. For x from the support Eag, the image fao{x]A) is contained 
in Eag by invariance of Eag. By continuity of fa, /a(x;A) C W for a close enough to uq. Since 
4>a depends continuously on a, 1 — a(o) < f)a{x)dm{x) converges to 0 as a —>■ oq. The 

proposition follows. □ 

6 Escape times 

In this section estimates for the average escape time from small neighborhoods of the support of a 
stationary measure that undergoes a bifurcation are derived, as function of the unfolding parameter. 

Let {{fa, da)}, a & I, he a, smooth one parameter family of random diffeomorphisms on M. 
Suppose that oq G / is a bifurcation value for {{fa, 9a)}- Let /x be a stationary measure for fag 
involoved in a bifurcation. Let VF be a small neighborhood of E. By Proposition 15.5L for W 
snfficiently close to E and a near ao, {{fa, ga)} possesses a nnique conditionally stationary measure 
fLa with support in W. Write for the density of fla- The transfer operator La acting on C^{W), 
has 4>a as a nnique fixed point. 

Let be the set of points in W with dVx^a H dW / 0 for x G X^. For i > 0, define 
_ f(^x^;dA). We suppress the dependence of W on a from the notation. 

Lemma 6.1 If Xi G X^, then for G f{xi‘,dA) one has Xi G i9I4i+i,a- 

Proof. This is clear from the definition. □ 

At X G X^, the boundary of Vx^a H W varies continuously but not smoothly with (x,a). It 
follows that La{(l)) cannot be expected to be more than continuous on X^ even for smooth (f>. 

Lemma 6.2 Suppose that (x,a) t—>■ f>a{x) is outside X^ U ••• \J X^~^, such that derivatives 
up to order k are bounded and their restrictions to a component of iy\(X'^ U • • • U X^~^) extend 
continuously to the boundary of the component. Then (x, a) La(j)a{x) is outside U • • • U 

X^. Likewise, derivatives up to order k + 1 are bounded and their restrictions to a component of 
IF\(X® U • • • U X^) extend continuously to the boundary of the component. 

Proof. For x 0 X^, the derivative of La(j) is of the form 

D{La(j)){x) = [ ■^ka{y,x)(l){y)dm{y) + [ na{y,x)ka{y,x)(l){y)dS{y). (22) 

jPc.anVF c'x d9(\4,„nM/) 

for a smooth function ka and a piecewise smooth function Ua (smooth outside the intersection of 
dW with dVx^a)- This identity shows that La4> is with bounded derivatives outside X^ , for any 
(j) G C^{W). The same holds for La{4>). Higher order derivatives are treated inductively. Similar 
to (ESI) one has 

d - f d f 

■^La(l){x)= ■^ka{y,x)(l){y)dm{y)+ Sa{y,x)ka{y,x)(t){y)dS{y) (23) 

-IV^c.anW Oa Jd{Va:,anW) 
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for some piecewise smooth function Sq. This shows that (x,a) La4>a{x) is outside for 
continuous functions {x,a) i—> (pa{x)- The derivatives of {x,a) i—> La(j)aix) on are bounded; 

moreover the derivatives on a component of extend continuously to the boundary of the 

component. 

Higher order derivatives are treated inductively. Suppose that (x,a) is outside 

U • • • U X^~^, such that derivatives up to order k are bounded and their restrictions to a 
component of U • • • U X^~^) extend continuously to the boundary of the component. Then 

(x, a) I—> La4>a{x) is outside X^ U • • - UX^. Likewise, derivatives up to order k + 1 are bounded 
and their restrictions to a component of 1T\(X^ U • • • U X^) extend continuously to the boundary 
of the component. 

The transfer operator La is the composition of the linear map La and the projection 


n(i?!)) = 4>/ [ 4'{x)dm{x). 

Jw 


The projection H is a smooth map which is well defined near cj)aQ in C^{W), a direct computation 
shows 

DU{(j))h = -— . h - 


Also, 


fw (f>(x)dm(x) {^^(t){x)dm{x)Y Jw 

5 * , k d' 


h{x)dm{x). 


/ (pa{x)dm{x) . , 

ocd Jw Jw 


■(j)a{x)dm{x). 


This implies that also (x, a) La{(t)a){x) is outside • -UX^ and has bounded derivatives. 

□ 


Proposition 6.3 For each k > 1, (f)a is outside • -UX^ ^ jointly in {x,a), the derivatives 

up to order k are uniformly bounded. 

Proof. Proposition 15.51 gives that fiaix) is continuous in (x,o). Recall that the support E of 
/i consists of finitely many, say k, connected components, permuted cyclically by the random 
diffeomorphism. An iterate of fag thus maps each component into itself. The restriction of the 
transfer operator to a small neighborhood of E has a single eigenvalue 1 and a remaining 
spectrum strictly inside the unit circle m- There is therefore no loss in generality to assume 
that the support of p. consists of a single connected component E, which we will assume for the 
remainder of the proof. 

Write = {Jj e C^{W) \ = Define the operator Ta : H°(W) H^{W) by 

Ta{lp) = Lai4>ao + "0) “ 4>ao- (24) 

Decompose (ja = 0ao + J’a, so that Ta{'4’a) = J’a- From the proof of Lemma HOI we get that Ta a 
smooth map on H^{W); 

DTa{fj)=DU{La{ct)ao+J^))La. 

However, La maps continuously differentiable functions to continuous functions, so that Ta does 
not define a map from H^{W), k > I, to itself. As a consequence we cannot obtain smoothness 
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properties of (pa by applying the implicit function theorem. To get smooth dependence of pa outside 
sets X® we reason as follows. We derive equations the derivatives of pa must satisfy, establish that 
the equations can be solved, and show that the solutions are the derivatives of pa- The reasoning 
follows the lines of the proof of Proposition IB.31 in Appendix I bI 

To prove that pa varies with a in points outside X^, note that -§^Pa{x) should be a solution 
Ma{x) to 

^Ta{Pa){x) + DTa{Pa)Ma{x) = Ma{x) (25) 

We claim that this equation is uniquely solvable. The spectral radius of DTa^ (0) = La^ is smaller 
than 1. As a consequence of the continuous dependence of La on a (see Proposition El, also 
DTa{pa) varies continuously with a. For a sufficiently close to oq, the spectral radius of DTa{pa) 
is therefore also smaller than 1. Hence 

OO 

(/ - DTa{Pa)y' =1 + Y,iDTa{Pa))\ (26) 

i=l 

see ES] This formula can be applied for DTappa) acting on C? functions. Indeed, DTa is compact 
on C C^{W), see Remark EH and has spectrum strictly inside the unit circle in C. From 

Ma{x) = {I - DTa{Pa))~^ ^TaiPa){x) 

= {l + DTa{Pa) + {DTaiPa)f + {DTa{Pa)f + ---)^Ta{Pa){x), (27) 

we get that Ma is continuous outside since it equals the sum of -^Tappa) and a uniform limit 
of continuous functions (compare Lemma |6.2I) . In particular Ma is uniformly bounded and has 
continuous extensions to the closure of components of VF\X^. We must show that 

\pa+h{x) - pa{x) - Ma{x)h\ = o{\h\) 

for X pL X^, as /i —?■ 0. Consider 'jaix) = pa+h{x) — pa{x) for x 0 X^. Then 

7a(x) = Ta+hiPa + iPix) - Taipa)ix) 

= DTa{paha{x) + ■^Ta{pa){x)h + R{x), (28) 

oa 

where 

R{x) = Ta+h{Pa +la){x) - Ta{pa){x) - DTa{pa)la{x) “ ^Ta(^a)(x)/l. 

We claim that for any e > 0 there is <5 > 0 so that \R\ < e(|7a| + \h\), if \h\ and |7a| are smaller than 
6. Since 7 a(x) is continuous in h we may further restrict 6 in this estimate so that |i?| < e(| 7 a| + \h\) 
holds for \h\ smaller than 6. Further, (/ — DTa{pa)) la{x) = ■^Ta{pa){x)h + R{x). Using (l2^ and 
the bound on |i?| gives | 7 a| < k\h\ for some k if \h\ < 6. Therefore |i7| < e{l + k)\h\ for some A; > 0, 
if \h\ < S. Now (1^ and (1^ give 

7a(x) - Ma{x)h = {I - DTa{pa))~^ R{x). 
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Using (P)|) it follows that \^a — Mah\ = o{\h\), /i —> 0. This proves that Ma equals the partial 
derivative ^^a- 

Higher orders of differentiability are proved by induction. Assume that (x, a) i—>■ 'ijja{x) has been 
shown to be outside U • • • U X^~^. Recall from Lemma Ih.21 that for maps (x, a) i—i- i/ja(x), 
■§^La{'4>a) is outside ■ -UX^. The right hand side of ( 1771 ) is therefore outside 
The above reasoning shows that Ma = -§^'4>a outside X^ U • • • U X^. Therefore -§^'4’a is outside 
X^ U • • • U X^. Also D-ipa = D{Ta{'ipa)) is outside X^ U • • • U X^, so that (x, a) i—> 'ipa{x) is 
outside X^ U • • • U The same clearly holds for (x, a) i—> ^a{x). □ 


Proof of Theorem ll.ldL We repeat the computation in the proof of Proposition 15.71 Let Ea 
be the interior of the support of ^a- Applying Lemma IHTH 


w 


Xa{x,uj)du°°{u)dm{x) > [ [ ^ dv°° {u)(j)a{x)dm{x) 

JEaJx* 4>a[x) 

> c [ [ Xa{x,ix>)dv'^{u)^a{x)dm{x) 

JEa JAfO 

c 

1 — a{a) ’ 


for (7 = 1/ max{x G W \ (j)a{x)}. By Proposition 16.31 (x, a) i—>■ ^a{x) is almost everywhere and 
has uniformly bounded derivatives. For each integer k there is a constant C with \4)a\ < C\a — ao|^ 
on W\EaQ. As in the proof of Proposition 15.71 we get that for each k there is a constant Ck > 0, so 
that 1 — a{a) < Ck\a — ao\^. □ 


7 Decay of correlations 

Consider a random family {fa} restricted to an isolating neigborhood IT of a stationary measure 
for all values of a form an interval I. The transfer operator La on Cq{W) possesses a single 
eigenvalue at 1. If the support of /Tq consists of r components, La has eigenvalues 0 < j < r, 

on the unit circle in the complex plane. These eigenvalues make up the peripheral spectrum of La, 
see Remark 12.51 In this section we consider bifurcations in which the number of components of the 
support of Ha changes. We will see how the rate of decay of correlations varies with the parameter 
a, providing a proof of Theorem ll.161 

Proposition 7.1 Let {fa}, I, he a family of random diffeomorphisms with an isolating neigh¬ 
borhood W. The eigenvalues and eigenvectors of the peripheral spectrum of La on Cq{W) vary 
smoothly with a. 

Proof. Let Xa be an eigenvalue that depends continuously on a and lies on the unit circle for 
a = oq. Since Hao is an isolated stationary measure, is a simple eigenvalue (see Remark 12.51) . 
Proposition iBRl in Appendix I rI implies the result. □ 
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Recall (HU and Lemma EH Write 

La<^ix) = / 

J A'* 

U^tpix) = / ij O fa{x-,U!i,. . . ,UJn)dl'{uJi) ■ ■ ■ dl^iiVn)- 

Ja” 

As in the computation for Lemma I2.1L 

/ Ul^p{x)'ip{x)dm{x) = / ip{x)U'^'ip{x)dm{x). 

Jm Jm 


( 29 ) 


After these preparations we now prove the statements on the speed of decay of correlations. First 
consider a single random map /. 

Proof of Proposition II . lH Let cj) be the stationary density. Write L'^tp = (/^ p{y)dm{y)) 4> + 
K^p. Compute 


tm[x] 


f p{x)U'^'ip{x)dm{x) = f L'^p{x)'ip{x)dr 
JM Jm 

p{y)dm{y) ) (/>(x) + 


'M 


'M 


'ilj{x)dm{x), 


so that 


'M 


p{x)U'^'4)[x)dm{x) — / p{x)dm{x) / '4)[x)cl){x)dm{x) 


IM 


<M 


IM 


p{x)il){x)dm{x] 


Note that the spectral radius of R is smaller than 1. By continuity of R, there is > 0 so that 
IIA^II < 1 for all a near oq. Hence for n € N, ||R’^|| < Crj"^ for some C > 0,r; < 1. The proposition 
follows from 


f R^p{x)'ip{x)dm{x) 

Jm 


< \\R'^'^\\c‘^{M)\\i’\\c'^{M) ^Crf'\\p\\c'i{^M)\\i’\\c'^{M)- 


□ 


Proof of Theorem czi This is proved by following the computation in the proof of Proposi¬ 
tion ^^31 above and noting that La has for a > no a single eigenvalue 1 and k — 1 eigenvalues that 
have moved smoothly into the unit circle. Write rja for the largest radius of the eigenvalues of La 
that lie inside the unit circle. As a consequence of the smooth dependence of the eigenvalues near 
the unit circle, see Proposition 17. II r]a is a smooth function of a. 

We claim that there exists C > 0 so that for all a near ao, ||R”|| < Crja- For a = ao, let E be 
the union of the eigenspaces for the eigenvalues in the peripheral spectrum. For a near oq, let Ea 
be the continuation of Ea^ = E. As Ea is finite dimensional and has a basis depending smoothly 
on a, it is clear that there exists C > 0 so that for a > oq, n G N, 

\\R^\Ej<Cri:. (30) 
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Let F be a subspace of C‘^{W) complementary to E. Write Pa for the projection to Ea along E. 
Then R = PaR + {I — Pa)R- By continuity of R and Pa, there is > 0 so that || ((/ — Pa)R)^ II < 1 
for all a near oq- Hence for n G N, 


{{I-Pa)RT\\ < 


(31) 


for some C > 0, < 1. Now pn|) and m prove the claim. As before, the proposition follows from 



R^ ip{x)'ilj{x)dm{x) 




□ 


8 One dimensional random maps 

The most complete description of bifurcations in smooth random maps is derived for random 
maps in one dimension. Consider a random endomorphism f{x',uj) on the circle S^. The random 
parameter oj is drawn from A = [—1,1]. What is proved below for random endomorphisms on the 
circle holds with obvious modifications for random endomorphisms on a compact interval that is 
mapped inside itself by all endomorphisms. 

A pathological example occurs if /(x; uj) is constant in x; the (unique) stationary measure is then 
a push forward of the measure on A. To avoid pathologies we assume the open and dense condition 
that the critical points of each map x i-H’ /(x;(j) have finite order. Also under this condition one 
finds that the regularity of stationary measures for random endomorphisms is substantially less 
then for random diffeomorphisms; their densities are only continuous. 

Theorem 8.1 The random endomorphism f G R°°{S^) possesses a finite number of ergodie sta¬ 
tionary measures //i,... ,//m with mutually disjoint supports Ei,... ,Em- All stationary measures 
are linear combinations of p-i,, Hm- 

The support Ei of Hi eonsists of the closure of a finite number of eonneeted open sets C^,... ,Cf 
that are moved cyclieally by /(•; A). The density fn of pLi is a funetion on 

Proof. The condition on the critical points of x /(x;u;) implies that 14 varies continuously 
with X. The reasoning used to prove Proposition 12.31 shows that the transfer operator L maps 
£^(S^) into C*^(S^). This implies continuity of invariant densities. □ 

Theorem 8.2 Let pL he an isolated ergodie stationary measure of f ^ R°°{§^) with density (f. 
Then eaeh f G i?°°(S^) sufficiently close to f possesses a unique ergodie stationary measure jl with 
support in V. The density of fi is close to cj). 

Proof. As in the proof of Theorem II.hi Note that L is a compact operator on C^(Al), compare 
the proofs of Theorems 12.31 and 15.31 □ 

Recall that iterates of a random map / are defined through ©• A periodic point x of pe¬ 
riod k is a point satisfying /^(x; wi,..., Wfc) = x for some uji,...,u>k G A. It is hyperbolic if 
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■^f^{x;(jOi,... ,Uk) at X = X differs from 0,1, —1. By the implicit function theorem, a family {/„} 
of random endomorphisms with fa^ = f possesses a hyperbolic periodic point Xa, Xao = x, for a 
near qq and for the same values ot iOi,... ,ujk, depending smoothly on a. 

Theorem 8.3 The set of stable random endomorphisms in is open and dense. 

Proof. Take / £ If the entire circle is the support of a stationary measure of /, then / is 

stable by Theorem l8.2[ Suppose that /r is a stationary measure whose support ill is a union 
of intervals Q mapped cyclically by /(■; A): /(Cj; A) = Cj+i (the indexes are taken modulo k). If 
pL is an isolated measure, / restricted to an isolating neighborhood of E is stable. 

The measure p is certainly isolated if for each boundary point x a E, either 

/^(x; A^) C interior E, 


or 

/^(x; cui,..., Wfc) = X, f^{x; iOi,... ,ujj) £ interior E 

for some uji,... ,uJk £ A, j < k. Indeed, invariance of E shows that in both cases /^(y; A^) £ E 
for any y near x. 

If not all boundary points are as above, then there is a boundary point x £ dE so that 
f{x] 0 Ji,... ,u)i) = X for I = k or I = 2k {I minimal) and {x;u!i,... ,ujj) £ dE for 0 < j < 1. 
Write xq = X and Xj = f^{x-,LJi,...,ujj) for j > 0. From Xj £ dE, Xj+i = f{xj;ujj+i) £ dE 
and 7^ 0, we see that Wj+i £ d/S.. Thus u}i,...,uji are all contained in 5A. Note that 

oji,..., uji) > 0 since otherwise x is an interior point of E. 

For / £ there are a neighborhood 17 of / and an integer N so that for each f ^ U, the 

support of the union of its stationary measures has at most N connected components. A random 
periodic orbit in the boundary of the support of a stationary measure of / £ C/ therefore has its 
period bounded by 2N. 

By transversality techniques a number of arbitrary small perturbations of / are carried through. 
The perturbations affect for cu £ dA and can be extended to other values of to using test 

functions. We will not present the detailed perturbations, but refer to m Section III.2] for a 
description of the techniques. By a small perturbation of / we may assume that the graph of each 
map /*(•; (5A)*), 1 < i < 2N, intersects the diagonal in x transversally. That is, 

(HI) the periodic orbits of period i < 2N for /(•; dA) are hyperbolic. 

There is then a bounded number of random periodic orbits with period bounded by 2N. A further 
small perturbation ensures that 

(H2) each periodic point x of period i < 2N is periodic for only one sequence 
(ji,..., Wj £ dA. 

Write V for the points in these periodic orbits. Recall that the number of critical points of /(•; dA) 
is hnite. A final small perturbation ensures that 

(H3) the critical values of /(•; dA) are disjunct from V. 
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Conditions |(H 1)| |(H 2)| |(H 3)| are clearly open and thus describe an open and dense subset of 


Consider / from this open and dense set. Let // be a stationary measure of / with support 
E. Let X be a boundary point of E belonging to a periodic orbit in dE. By[(in)| X belongs to a 
hyperbolic periodic orbit. By (H 2)[ there is a unique graph ... , 0 Ji) with wi, ■ ■ ■ , 0 Ji G dA 

through X. It is not possible that~^/^(x; wi,..., toi) > 1, since other orbits would then be repelled 
and X would not be in the boundary of E. Hence 0 < ... ,uJi) < 1: the random periodic 

orbit through x is an attracting periodic orbit for ... ,uJi). By |(H3)| there are no interior 

points in E being mapped onto x under iterates of /. As a consequence, /i is isolated. Therefore / 
is stable. □ 


As a next step we consider one parameter families of random maps and show that bifurcations 
typically occur at isolated parameter families. The theorem below moreover describes the possible 
codimension one bifurcations. The space of smooth families of smooth random maps x fa{x\uj), 
a ^ I, will be given the uniform topology as maps {x,oj,a) i—> fa{x\u)) on x A x /. 

We start with a description of three types of bifurcations caused by violation of one of the 
conditions |(H 1)[ |(H 2)[ [(H 3)| These are proved to be the only codimension one bifurcations. 

Random saddle node bifurcation Random homoclinic bifurcation 




Figure 2: Consider a random map fa^ for which points are mapped randomly into the region bounded by 
the two graphs. Depicted on the left are the graphs a random map /(you), uj G 9A, with a random saddle 
node bifurcation. The support of the stationary density is the interval between the hyperbolic fixed point 
of the lower map and the saddle node fixed point of the upper map. The right picture shows graphs of a 
random map with a random homoclinic bifurcation. Here the support of the stationary density stretches 
from the left hyperbolic hxed point of the lower map to the critical value of the upper map. 
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Definition 8.4 The smooth one parameter family of random endomorphisms fa on the cirele un¬ 
dergoes a random saddle node bifurcation at a = oq, if there exists x in the boundary of the support 
of a stationary measure such that 

= X, -^fa^{x]UJi,...,UJk) = l (32) 

for some ui,... ,u}k £ dA. The random saddle node bifureation is said to unfold generically, if 

(^) ^0’ ■^faix-,(jJi,...,uJk) 7^0 (33) 

at a = oq. 

Definition 8.5 The smooth one parameter family of random endomorphisms fa on the circle un¬ 
dergoes a random homoclinic bifurcation at a = oq, if there exists 

• a stationary measure fi with support E with a hyperbolic periodic point Xa in the boundary of 
E for all a near oq, and 

• a critieal point ya for fa{-',uji), oji G dA, in the interior of E, 
sueh that 

/ao(yao;wi,...,u;0 = Sao (34) 

for some u) 2 ,..., 0 Ji £ dA. The random homoelinie bifureation unfolds generically if 

^ (/a(ya;wi,. . . - Sa) 7^ 0 (35) 

at a = qq. 


Definition 8.6 The smooth one parameter family of random endomorphisms fa on the cirele un¬ 
dergoes a random boundary bifurcation at a = oq, if there exists x in the boundary of the support 
of a stationary measure and (wi, ... ,uJk) 7 ^ (wi, ... ,0^) £ {dA)^, such that 

= X, -^fl;^{x;uJi,...,uJk) e {0,1) (36) 

and ^ 

faoi^',‘^i,---,i^k) = x, —/^Q(x;a;i,... jWfc) £ ( 1 , 00 ) (37) 

Write Xai^^i, ■ ■ ■, ^k) and Xa{Cbi,..., Ok) for the eontinuations of the hyperbolic periodic points. The 
random boundary bifurcation is said to unfold generically, if 


d d 

{Xa{a^l: • • • j OJk)] . . . , LOk) 7^ iXa{Ol, ■ ■ ■ , Ok] )Oi, . . . , Ok) 


( 38 ) 


at (X — qq . 
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Random boundary bifurcation 



Figure 3: Depicted are parts of the graphs of a random map f{-;uj), uj £ dA. The solid curves lie on two 
of the graphs of /(/(•; uji); 1 ^ 2 ), ^ 1 , 0^2 £ dA, intersecting in a point that lies on two hyperbolic period two 
orbits (one stable, one unstable) distinguished by different uj values. A random boundary bifurcation results 
if this point lies on the boundary of the support of a stationary measure. 


For an open interval I, write R^{I,M.) for the space of families of random maps in R^{A4) 
depending on a parameter in I. Equip the space R^{I,A4) with the topology. 

Theorem 8.7 For {fa} from an open and dense subset 0 /S^), fa has only finitely many 
bifureations. A bifureation point is a random saddle node bifurcation, a random homoclinic bifur¬ 
cation, or a random boundary bifurcation and is generically unfolding. If the number of stationary 
measures is loeally constant at a bifurcation point, the bifurcation is an intermittency bifurcation. 
Otherwise the bifurcation is a transient bifureation. 

Proof. For a bifurcation value for a family in S^), either [(H 1)[ |(H 2^ or ME is violated. 

Similar transversality arguments as in the proof of Theorem 18.31 show the following. For an 
open and dense subset of R°°{I,§^), at most one of these conditions is violated at a bifurcation 
value and the resulting bifurcation unfolds generically as stated in Definition 18.4118.hi or 18.til Since 
the random bifurcations are unfolding generically, they occur isolated. □ 

9 Case studies 

In this section we illustrate the general theory on two examples; a randomized version of standard 
circle diffeomorphisms and a randomized version of logistic maps on the interval. We explain how 
random saddle node bifurcations occur in both examples and random homoclinic bifurcations in 
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random logistic maps. For the random circle difFeomorphisms we consider rotation numbers and 
study their dependence on parameters. The reader is referred to m for the theory of deterministic 
circle and interval maps. 


9.1 Random circle diffeomorphisms 

The standard circle map acting on x € M/Z and depending on parameters a,e is given by 

fa{x) = X + a + — sin(27rx) mod 1. 

ZTT 


Consider fa for a fixed value of e G (0,1) for which fa is a diffeomorphism. Introduce the lift 
Fq : M I—> M, 

Fa{x) = X + a H-sin(27rx). 

27r 

It is well known that the rotation number pa of fa, 


Pa = 


lim 

k—^oo 


F^{x) — X 

k 


(39) 


is well dehned and independent of x. The rotation number depends continuously on a. The 
rotation number is rational precisely if fa possesses periodic orbits. For a hxed rational number r, 
the rotation number of fa equals r for an interval of a values. In the interior of such an interval, 
fa has exactly one hyperbolic periodic attractor and one hyperbolic periodic repeller, see ESI. 

In the following we consider standard circle diffeomorphisms with a random parameter: 

/a(x; u;) = X + — sin(27rx) + a + aco (40) 

27r 


for X G M/Z and a random parameter oj chosen from a uniform distribution on A = [—1,1]. The 
value of a determines the amplitude of the noise, we assume it has a hxed value. We consider hxed 
e G (0,1) for which x i—> fa{x;Lv) is a diffeomorphism. Write 


Fa{x‘,uj) = X + a + aoj + — sin(27rx) 


(41) 


for the lift of /a(x;a;). Note that Fa(x;(j) — x is periodic in x with period one. 


Proposition 9.1 For each parameter value a, the random standard circle family {fa} has a unique 
stationary measure pa- The density (j)a of pa is smooth and depends smoothly on a. The support 
of Pa is either the entire eircle or finitely many intervals strictly contained in the circle. The latter 
possibility is only possible if pb is rational for eachb G [a — a,a + a]. Bifurcations where the support 
of Pa changes discontinuously, are generic saddle node bifurcations. There are finitely many such 
bifurcations. 


Remark 9.2 Observe that fa has a hyperbolic fixed point for a G (“^)^)- Hence, fa has a 
stationary measure supported on a single interval precisely if both o — a > and a + a < ^. 

This occurs for a nonempty interval of a values if cr < -f-. 
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Proof. It is well known that a circle diffeomorphism with irrational rotation number has its 
orbits lying dense in M/Z. It follows that if the family of circle maps fa{x;Lv) for varying u € A 
contains a member with irrational rotation number, there is a (necessarily unique) stationary 
measure supported on all of M/Z. 

Suppose now that fa{x] u>) for each uj € A has rational rotation number Pa+au = p/q- Write x^ 
for a periodic point from a periodic attractor of x fa{x-,uj) depending continuously on tu. Recall 
that X I—> fa{x;io) has a unique periodic attractor. Let I 4 = yj^i^AXui = [x_i,xi]. The random 
standard family is increasing in x and in w, so that for all x G 14 ; and all u G we have 

x-i = i4(x_i;-l) < F^{x-i;u) < F^{x;u) < T„^(xi;a;) < F«(xi;l) = xi. 

It follow that the orbit of 14 is invariant. For a fixed w G A, all points outside the unique periodic 
repeller of fa{-;ix) are attracted to its periodic attractor. This implies that there is a unique 
stationary measure supported on the orbit of I 4 . 

Compute 

■^faiF ^)=Y1 bfuj] u)). 

a oa X 

As all terms in the sum are positive, a random saddle node bifurcation occurs isolated. The random 
family {fa} therefore has only a finite number of random saddle node bifurcations. □ 



Figure 4: Numerically computed stationary densities of the random standard circle map. On the left for 
\a\+a < e/27r, on the right for |a| + cr > e/27r. The explosion of the support of the stationary density follows 
a random saddle node bifurcation. 

We define the rotation number for the random standard circle map, when its exist, by 

(42) 

k—*oo k 

The rotation number measures the average rotation per iterate of fa- Note that pa is a random 
variable, depending also on the starting point x. 
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A simple but usefull lemma shows that pa is independent of the initial condition x. 


Lemma 9.3 If pa{x;u>) exists for some x G M/Z,a; G then x i—> pa{x\u> exists for all x G M/Z 
and is constant in x. 

Proof. Observe that is a lift of /^(•;cj), so that F^{x](jj) — x is periodic in x with period 

1. Thus 

maxiFN x: io) — x| — min{F„^(x;a;) — x} < 1. 

sGR ^ ^ a;gR L a \ I / 

Compute 

\F^{x-,u:)-F^{y;u)\ < \{F^{x-,u) - x) - iF^{y,uj) - y)\ + \x - y\ 

<l + \x-y\, 


so that 


lim 

fc—>-O0 


Fa{x]^)-x F^{y;u)-y 


k 


k 


= 0 . 


It follows that the limit linifc^oo ^ ; if h exists, is independent of x. 


□ 


Write 


Fa{x]u) = X + (5a(x;u;), 


where the function da(x;Lo) is periodic with period one in the variable x. We can consider <5 as a 
function defined on M/Z. A simple induction argument gives for each A: G N, 


k-l 


faix;uj) = X + ^(5o s"{x-,u) 


(43) 


i=0 


where S is the skew product system (see equation @) on M/Z X . Recall that pa x 1^°° is an 
S'-invariant measure. 


Proposition 9.4 


Pa — 


K{6{x]uj))dpais) — a.s. 


(44) 


where E is the expectation operator. The right hand side of is independent ofx and is a smooth 
and nondecreasing function of a. 


Proof. Consider the Birkhoff sum in equation (ESI) 

f^{x-,uj)-x 1 


k 


o S''{s-,Uj). 


i=0 
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Devil’s staircase {a = 0.00) Smooth devil’s staircase (cr = 0.05) 



Figure 5; The function a ^ pa- On the left the devil’s staircase; the rotation number of the deterministic 
standard family. On the right the rotation number of the random standard family. 


By Birkhoff’s ergodic theorem, 

- X 


lim 

k^oo 


k 


/ / f5a{s-,u:)p,a{ds)v°°{du;) p,a 

J J 

/ E(/3a(s;a;))/ra(ds) Ha x 1 ^°° - a.s. 

Jr/I, 


(45) 


The fact that the rotation number "when it exist is independent of the initial point, see (Foil . implies 
that this equality holds for all x and v°° — a.s.. 

Write a i—> h{a) for the right hand side of (I44j) . Smoothness of h follows from smoothness of 
the stationary density (pa and (1^ . Write the standard family as Ra o /(.;a;) where is the 

random map f{x;Lo) = x + sin(27ra;) + cr^{Lo) and Ra is the translation with coefficient a. Then 
for ai < 02 and k > 1, {Rai ° < {Ra 2 ° and thus, 


Pai = lim 

fc—>cx: 


{Ra,0f{.-U)f-Id 

k 


< lim 


{Ra2 0f{.-U)f -Id 


k 


Pa2 ■ 


□ 


9.2 Random unimodal maps 

This section is devoted to the investigation of the randomized version of the logistic family 

fa{x;uj) = {a +aLo)x{l-x). (46) 

on [0,1]. The random parameter uj be chosen from a uniform distribution on A = [—1,1]. Through¬ 
out this section we will assume that 

fl + crcj G (1,4), (d7) 
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for all w E A. As a consequence, the interval [0,1] is mapped into itself by each map x i—> fa{x;u;) 
and the fixed point at the origin is repelling. Any stationary measure will therefore have support 
contained in (0,1). We will demonstrate that there is only one stationary measure. 

Proposition 9.5 The random logistic map fa has a unique stationary measure. 

Proof. We collect some facts from unimodal dynamics needed in the sequel of the proof. The 
following facts hold for unimodal maps with negative Schwarzian derivative such as the logistic map 
X I—> ax{l — x). By Guckenheimer’s theorem, see |321 Theorem III.4.1], x ax{l — x) possesses 
a unique attractor Aq. The attractor Aq is either a periodic attractor, a solenoidal attractor, or 
a finite union of intervals on which the map acts transitively. In all cases, the omega-limit set 
of the critical point c (with c = | for the logistic map) is contained in A^. In fact, if A^ is not 
a periodic attractor, then c is contained in Aq. It follows from a result of Misiurewicz, see m 
Theorem III.3.2], that the basin of attraction of Aq is an open and dense subset of (0,1). 

We will distinguish the following two cases. 

Case (i): There exists cu E A, so that c E ^a+auj-, 

Case (ii): otherwise. 

The two cases are treated separately. 

Case (i): Write 

W=f][Jf[{^) 

n>0 i>n 

for the omega-limit set of c under all possible random iterations. Observe that W is an invariant 
set. From the properties of the noise, W consists of a hnite union of intervals. Note that c E W, 
so that W equals the closure of the positive orbit Ui>o/a(c;^^) of c under all possible random 
iterations. We will prove that for each x E (0,1), y E W and e > 0, there exist n > 0 and u E A^ 
with the property that 

-y\ < £■ (48) 

This implies the W is the unique minimal invariant set, which in turn implies the theorem in the 
first case. 

Fix X E ( 0 , 1), y E W, e > 0 . From the construction of IF, there exist cji E A^, i > 0 , so that 
|/*(c; o^i) — y| < e. By continuity of x i—> /*(x; o^i), the same holds with c replaced by a point from 
a 5 neighborhood of c for some d > 0 . We need to establish the existence of A E A^ and j > 0 so 
that I/a (x; A) — c| < 5 . Let u;2 E A be such that c E ka+(Tui2 ■ Since the basin of attraction of Aa+o-ojj 
is open and dense, there exists cua E A with xi = fa{x] 0 jf) contained in the basin of attraction of 
Aa+(Tta2- For i large, Xj = /*“^(xi;a;2,x’2, • • ■ ) is as close as desired to ka+cruj2- If ^a+(juj2 i® ^ finite 
union of intervals, we get that Xj is contained in Ka+(juj2 for large enough i. As inverse images of 
c for X I—> fa{x](jJ2) are dense in Aa+auj2^ one deduces that there exist U4 E A^ and k > 0 so that 
/^(xi+i;a^4) lies in a 5 neighborhood of c. Indeed, if Aa+auj2 i® a finite union of intervals, then we 
find W5 with Xj+i = f{xi;u}^) equal to an inverse image of c. If Aa+cruj2 i® a solenoidal attractor, 
then X I—>■ fa{x\uj2) is infinitely renormalizable. In this case one can use 0^4 = {u>2,a)2, • • • )• Also 
for a periodic attractor Aa+o-ta2 containing c one uses 0^4 = (u;2,to’2, • • ■ )• 
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Case (ii): By Guckenheimer’s theorem, x fa{x;io) possesses a unique periodic attractor for each 
a; G A. Write V for the the union of Aa+au over lo € A. Dehne 

i>0 

This is clearly an invariant set. Note that we do not claim that c is outside of W. Arguments as 
before prove with this definition of W: for suitable noise one finds an orbit starting at x that 
approaches a point in V and then with further iterates approaches y ^W. □ 

Remark 9.6 The above proof applies to show that a random unimodal map g{x;u;) with negative 
Schwarzian derivative for each u) (the invoked theorem by Guckenheimer is true for these maps) 
has a unique stationary measure. 



Figure 6: Numerically computed Birkhoff averages lim„^oo of fix) = l[o. 4 , 0 . 6 ] for the 

logistic family (left picture) and the random logistic family with a = 0.005 (right picture), for parameters a 
ranging from 3.8 to 3.9. The flat part in the left picture, where the average equals 1/3, runs from a saddle 
node bifurcation to a homoclinic bifurcation. The numerical computations show that for the random logistic 
family these are replaced by their random versions; for parameters from the flat part the stationary measure 
is supported on three disjoint intervals cycled by the random map. 


Perturbing away from the deterministic logistic family one sees that both random saddle node 
bifurcations and random homoclinic bifurcations occur in the random logistic family {fa} for small 
noise levels. Typically one can expect the following scenario. We start by recalling some facts 
concerning the dynamics of the deterministic map /a(-; 0). The map /a(-; 0) is called renormalizable 
if there exists an interval I and a positive integer q, so that fail] 0) C I. Let [a_.a+] be a maximal 
interval so that /a(-;0) is renormalizable for a G [a_,a+] with q constant. Then fa undergoes a 
saddle node bifurcation at o = a_ involving a periodic orbit of period q. At a = o+, fa undergoes 
a homoclinic bifurcation, where an iterate of fa maps the critical point onto a periodic orbit of 
period q. For small noise levels (i.e. a small) one expects a random saddle node bifurcation near 


37 










a = a- and a random homoclinic bifurcation near a = a+. Figure IHl illustrates this by computing 
Birkhoff averages for the logistic family and the random logistic family. See m for explanations 
of the computations for the logistic family. 

A Representations of discrete Markov processes 

In this appendix we explore the relation between random maps and discrete Markov processes given 
by stochastic transition functions. The random maps considered in this paper depend on random 
parameters, where the number of random parameters equals the dimension n of the state space Ai. 
Proposition lA. II gives a wide class of Markov processes that can be represented by random maps by 
n random parameters. The Markov process given by random maps depending on a larger number 
of random parameters (or even given by some measure on the space of maps) can be represented 
by random maps with n random parameters. 

Iterating a random map involves more random parameters obtained by independent draws at 
each iterate. By means of an example we explain how random maps with a smaller number of 
random parameters may be brought into the context of this paper. Consider the delayed logistic 
map Xn+i = — Xn-i)- Let Vn+i = Xn- This defines a dynamical system {xn+i,yn+i) = 

{yxn{^ — yn),Xn)- Assume now that /r is a random parameter varying in some interval with some 
distribution. This yields a random map 

f{x,y;i^) = {yx{l - y),x). 

The derivative Df is singular along x = 0. As // is a single random parameter, this random 
diffeomorphism does not fit into the context considered in this paper. Considering two iterates 
gives two independent draws (/z, z/) of the random parameter (that is, random parameters taken 
from a square) and yields the random map 

f{x, y, fj,, v) = {y,x{l - y), vyx{l - x){l - y)). 

If X and y stay away from 0 and 1, the map and the dependence of (/z, v) are injective. The second 
iterate of the delayed logistic map with bounded parametric noise fulfills the assumptions in this 
paper. 

There are other examples of maps with parametric noise that cannot be made to fulfill the 
assumptions used in this paper. For instance, random maps f{x;uj) = x + {x — w)^ with random 
CO from an interval, fail to satisfy the injectivity assumption of cj i—;• f[x]co). If oo is chosen from 
a uniform distribution, then the density of the transition function will not be bounded. Figure [3 
indicates a random boundary bifurcation for a similar random map. 

Consider discrete Markov processes given by transition functions P{x,-). The following prop¬ 
erties hold. 

For fixed A £ B, x T’(x, A) is measurable. 

For fixed x £ M”, P{x, ■) is a probability measure. 
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Random boundary bifurcation 



Figure 7: A random map f(x\w) = f{x — w; 0) + w with a random boundary bifurcation of the stationary 
measure with support between the ordinates indicated by dotted lines. If u is chosen from a uniform 
distribution, then the density of the transition function will not be bounded. 

Denote by y i—> k(x,y) the density of P{x,-). Write Ux for the support of k{x,-) and let U = 
yjx ({rr} X Ux)- Assume that Ux is diffeomorphic to the closed unit ball A in M” and varies smoothly 
with X. We will assume that y k{x, y) depends smoothly on (x, y) £ U, meaning that k can be 
extended to a smooth function defined on an open neighborhood of U. Under these conditions we 
will construct a representation by a finitely parameterized family of endomorphisms. That is, we 
will construct a family of endomorphisms {/^} on M"", with parameters /i from an n dimensional 
ball, and a measure v on the parameter space so that P{x,A) equals £ A | /^(x) £ A}. A 
corresponding result holds for discrete Markov processes with noise from an n-dimensional box, see 
m Appendix D]. See [30] for a discussion of the existence of representations by sets of measurable 
or continuous maps. The paper im contains a result on representations by endomorphisms, under 
the assumption of unbounded noise. 

Proposition A.l There is a family of endomorphisms f^, /i £ A, and a measure n on A with 
smooth strietly positive density, so that 

1. {x,y) f^i{x) is smooth, 

2. for each x £ At, // i—> /^(x) is injeetive, 

3. P{x,A) = v{y^A \ /^(x) £ A). 

Proof. We follow the arguments in HD, combined with the use of polar coordinates to map the 
unit ball A to [0,1]”. Let : 14 —> A be a diffeomorphism, depending smoothly on x, from the 


39 






support Vx oi y 
ball, 


k{x, y) to the unit ball. Consider polar coordinates x ^ [0,1]"' 


x(6) • • • 5?n) = 6 


cos(7r^2) 

sin(7r^2) cos(7r^3) 

sin(7r^2) • • • sin(7r^„_i) cos(27r^„) 

\ sin(7r,f2)---sin«„_i)sin(27r^„) ) 


A on the unit 


For ^ G [0,1]"', define sets 0 < z < n, by 


i n 

j=l j=i+l 

Write Ci{C) = V’^ and let lo = (wi, ...,ujn) be given by 


Ui = 


'CiiO 


k{x,y)dm{y) 



Ci-iiO 


k{x,y)dm{y) 


Since A: > 0, a; = ©(^ gives a 1-1 correspondence. Let ry = J^,^^~^dm{y) j ^^^~^dm{y). Here 
r/ = is a 1-1 correspondence. The correspondence a; —> t/ is a smooth diffeomorphism as k is 
smooth and strictly positive. Then 


ff,{x) = ^©X \f^)- 


gives the required smooth random maps. 


□ 


For discrete Markov processes on a circle there is an easy necessary and sufficient condition on 
the transition maps for a representation by random diffeomorphisms. 

Proposition A.2 Let A4 he the circle endowed with Lebesgue measure. Write 14 = [^-{x),l+ix)]. 
There is a representation by random smooth diffeomorphisms if and only if 

Q 

-k{x,l-{x))l'_{x) + I —k{x,y)dy^0 

Ji-{x) dx 

for z G 14 . 


Proof. The construction of the representation by random smooth maps proceeds as follows. For 
^ G [0,1], write C'(^) = [l-{x),l-{x) -|-^(Z+(x) — Z_(x))] and let 


Since A: > 0, the map ©, ©(^ 
phisms is given through 


OJ = 


k{x,y)dy. 


Jc{0 

OJ, is a diffeomorphism. The representation by random diffeomor- 


f^{x) = l-{x) + G ^{oj){l+{x) - l-{x)). 
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Note that for fixed u 


jL 

dx 


k{x,y)dy = 0. 


(49) 


J[l_(x),fU^)] 

With, say, I- < (IHHl yields 

-k{x,l-{x))l'_{x)+ k{xjuj{x))fi{x) + [ -^k{x,y)dy = 0. 

Ji-{x) ox 

Hence f^{x) = 0 precisely if -k{x,l-{x))l'_{x) + ■§^k{x,y)dy = 0. 


□ 


B Regularity of solutions of integral equations 

In a number of places in this paper eigenvalue equations Lacj) = Ac/) for the transfer operator La 
arise. As La depends only on a, a direct application of the implicit function theorem as found 
e.g. in nniiis] yields only weak regularity properties of the solutions. 

The following remark is a variant of Proposition 14.11 It allows an application of m Proposi¬ 
tion 3.6.1] to show that eigenvectors and eigenvalues of La, in the case of simple eigenvalues, vary 
smoothly with a. 

Remark B.l For a £ I, a La is a map from I into the space of 

bounded linear maps from into C^{A4). 

We include in this appendix an alternative route to obtain such smoothness, as introduction to 
the more involved reasoning in Section 121 

Given is Laofao = ^ao4>ao with La acting for o near oq on C^{A4) (here we are considering 
complex valued functions). Similarly we can consider La acting on Cq{W) for an isolating neigh¬ 
borhood W. Assume that Aqq is an isolated eigenvalue of La^- Denote by E the span of fao and 
let F^(IT) be a complement of E in C^{M). Consider functions fia = 4>ao E'f’a with ifa G F^{M). 
We wish to solve Lafa = Let P be the projection to E along F^{A4). Considering a second 

parameter A, La4>a = Ac/)^ decomposes as 

f {I - P)Lai4>ao +i’a) = AV’a, 

\ PPaifao T V’a) — A(/)(iq . 

The top equation can be solved for fja as a function of A and a for a near oq and A near Aao- 
In fact, by the Fredholm alternative, ipa = ((L — P)La — A/)”^ ((/ — P)Lafao)- Putting this into 
the bottom equation yields a single equation for A. Note that for stationary measures, A = 1 
automatically solves this equation, compare the proof of Theorem 11.101 in Section ^ 

Write the top equation as a fixed point equation Ta'il’a = 'f’a with parameters a. The map Tq, 
is a compact linear map mapping F^{M) into compare the proof of Proposition 12.3L 

Lemma B.2 If (x,a} i-^ ipalx) is then {x,a) ^ Taifaix) is 

Proof. See Section El □ 
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Proposition B.3 Consider the integral equation T^ipa = V’o with Tq, as above. The fixed point 
X I—> fia{x) is smooth jointly in x,a. 

Proof. Given is a unique fixed point ipa depending continuously on a. Formally differentiating 
Tafi>a = fi’a witli respect to a gives 

d d d d 

^ (Tafiaix)) = -^Tafia{x) +T„ —= —fia{x). 

So should be the solution Ma of -^Taipa{x) +TaMa{x) = Ma{x). That is, 

= (50) 

By the Fredholm alternative m, I — Ta is invertible on F^{A4). The right hand side of (1501) 
is therefore a continuous function. To establish that is the derivative we must show 

\ijja+hix) —fia{x) — Ma{x)h\ = o(|/i|) as /i —> 0 (compare the proof of the implicit function theorem 
in e.g. m or [El)- Write 7 a(x) = fia+h{x) — 'tpa{x). Now 

7 q ,( x ) = Ta+hilpa + 7a){x) - Tai'4!a){x) 

= Tajaix) + ■^Tafia{x)h + R{x), (51) 

oa 

where R{x) = Ta+hifila + 7a)ix) - Ta{'lpa)ix) - Taja{x) - -§^Ta1pa{x)h. Since Ta{fi) is 

differentiable, for any e > 0 there is 5 > 0 with \R\ < e(| 7 o| + \h\) if | 7 q,|,|/i| < 6. Since 7 ^ is 
continuous in h, we may further restrict 5 so that this estimate on |i?| holds for \h\ <6. From (j51 () 
we get 'Ja = {I — Ta)~^{-^Ta'i’ah + R) so that | 7 q,| < C\h\ for some C, if \h\ < 6. This implies 
\R\ < e(l + C)|h| for \h\ < 6. As 


(/ - rQ,)(7„ - Mah) = R 

(from (ISni) we get (/ — Ta)Ma = -^Tafia), we derive | 7 o(x) — Maix)h\ < Ke\h\ for some AT > 0, if 
\h\ < 6. This proves that equals the partial derivative 

Higher order derivatives are treated by induction. Assume that {x, a) ipaix) has been shown 
to be Cfi By Lemma FB.21 {x,a) 1 —> Tafi>a{x) is So Dfia = D{Ta'fia) is Cfi As (/ — Ta)~^ 

maps to F^(J\A), the right hand side of (15011 is a function. The above reasoning shows 

that Mq, = ■§^'Pa- Therefore ■§^'4>a is Cfi so that (x,a) 1 —>■ fiaix) is □ 
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